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METHOD AND DEVICE ARRANGEMENT FOR AUTOMATIC DOSE 
CONTROL OF CHEMICALS 

5 

FIELD OF THE INVENTION 

The present invention relates to a method for automatic dose control of chemicals in 
liquid treatment system, particularly for automatic dose control of chemicals in 
10 water, waste water and industrial water treatment. 

The present invention also relates to a device arrangement for automatic dose 
control of chemicals in liquid treatment system. 

15 

BACKGROUND OF THE INVENTION 

Several short and long-term phenomena appearing in liquid treatment processes 
influences dosing of liquid treatment chemicals. There are often large variations in 
20 the quantity of liquid to be treated. For example the consumption of industrial or 
drinking water and the amount of waste water greatly vary depending on the plant 
capacity, time, season, weather, etc. The quality of the incoming water, e.g. 
turbidity, suspended solids, pH, phosphorus, temperature, also varies rapidly 
sometimes. 

25 

Today, chemicals are often dosed on the basis of incoming flow rate in liquid 
treatment plants. This does not guarantee the treatment efficiency, if the quality of 
liquid greatly changes. This leads easily to over or under dose of chemicals 
depending on the quality of liquid. Many complex phenomena, such as coagulation 
30 and flocculation take place in liquid treatments. The effect of chemicals on 
treatment result is non-linear. 

Feedforward dosing on the basis of incoming flow rate does not ensure sufficient 
treatment results since the quality of the liquid is changing. Feedback control with a 
35 PID (Proportional, Integral and Derivative) controller provides improvement if 
fluctuations are fairly small. Variations in operating conditions require adaptation of 
the controllers. A gradual adaptation of the PID parameters can be done by gain 
scheduling but it is limited to a fixed control structure. Flexible control structure can 
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be realized by switching strategies. However, both gain scheduling and switching 
strategies require complex logic for adaptation in multivariable systems. Also 
smooth and reliable operation is difficult to achieve. The PID parameters can also 
be tuned on-line, e.g. by self-tuning PID controllers. 

5 

One example of adaptively designing self-tuning PID controller is disclosed in 
WO0198845 and it can be characterized by parameter values derived from 
interpolation of process model parameters. Parameters characterize each of the 
models. A value of the parameters is selected from a set of predetermined 

10 initialization values. For each parameter value so-called accumulated Norm (which 
is derived from a model square error calculated for the models) is calculated as 
repeated evaluations of models are conducted. For each parameter an adaptive 
parameter value (a weighted average of the initialization values) is also calculated. 
The set of adaptive process parameter value are then used to redesign a process 

15 controller. Parameterized controller is applied to adaptive feedforward/feedback 
PID control. 

Self-tuning PID controllers are efficient in adapting only to slow changes in 
operating conditions, as such adaptation mechanism is fairly slow. Disturbances 
20 during the adaptation may cause serious problems. Limitations become more 
prominent in multivariable systems. 

Adaptive fuzzy control is less sensitive to disturbances than self-tuning PID 
controller. One example of multi-region fuzzy logic control system is disclosed in 

25 US 6041320. The system utilizes an auxiliary process variable to determine which 
of several regions of different gain a non-linear process is operating. The operation 
area of a process is divided into different regions in advance. Each region has its 
own controller which is implemented by a linear set of membership functions and 
rules. The auxiliary variables in each rules is used for selection of proper region. 

30 Solution becomes impractical with increasing number of regions and auxiliary 
process variables as each region has its own controller. This leads large number of 
parameters and difficulties to generalize the control system to new applications. 

All these methodologies lead to extremely complicated control structures containing 
35 multiple interactive and dependent parameters which are difficult to handle. The 
adaptation logic becomes complicated as the number of parameters handled 
increases. The operation of different controllers can be presented with nearly similar 
control surfaces in certain process conditions. However, it is essential how the 
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adaptation mechanism is done. Compact adaptation structure is necessary in 
changing operating conditions and implementing controllers into new processes. 
The present inventors have published a conference paper concerning general state of 
the art (Esko Juuso, Katja Virret and Marjatta Piironen "Intelligent methods in 
5 dosing control of water treatment", in: Proceedings of Workshop on Applications in 
Chemical and Biochemical Industry, ERUDIT, September 15, 1999). Paper 
describes a preliminary modeling experiment for water treatment by linguistic 
equation (LE) method but not LE-based control in water treatment. In the paper it is 
mentioned that LE-controller can be adaptive like any other controller. 

10 

The current invention is founded on a novel approach where dosing control is based 
on compact control system with adaptive and non-linear features. Smooth and fast 
adaptation to rapid variations improves the performance of dosing control. Adaptive 
LE controllers can operate at large range of process conditions because the 

15 adaptation is based on detection of the process state in advance. Predefined 
adaptation approaches addressed by the inventive method for dose control does not 
require time consuming identification of models or parameters. Compact parameter 
set is also beneficial in the tuning of the LE controllers. In the present invention the 
number of parameters is small compared to the current adaptive controllers known 

20 in the art, and all the parameters are understandable giving insight on the process 
operation. 



SUMMARY OF THE INVENTION 

25 

The present invention provides a method and a device or device arrangement for 
dose control of chemicals. This is accomplished by using the properties of liquid to 
modify the control surface of a linguistic equation (LE) controller adaptively by 
means of a predefined adaptation model, e.g. a model developed before the use of 

30 the control, to control the dosing of one or more chemicals to the liquid by one or 
more said controllers. The linguistic equation (LE) controller already handles non- 
linear operating conditions. Non-linear features are modified with non-linear 
transformations. Adaptation further extends the operating range of controllers 
without changing the control structure which comprises LE controllers and 

35 adaptation models. Therefore adaptive LE controllers can operate at large range of 
process conditions. 
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It is characteristic for the invention what is claimed in the independent claims. 
Some of the preferred embodiments are claimed in the dependent claims. 

It is an object of the present invention to provide a new method for dose control of 
5 chemicals. 

It is another object of the present invention to provide a new device and device 
arrangement for dose control of chemicals. 

10 It is still another object of the present invention to provide a new method for dose 
control of chemicals in situations where the environment, such as the quality or the 
flow rate of the liquid, is changing rapidly. 

It is still another object of the present invention to provide a new method for dose 
1 5 control of chemicals in water treatment. 

In the present invention the amount of one or more chemicals are monitored and 
controlled with non-linear adaptive multivariable controllers on the basis of 
measurements in liquid treatment system. Preferably the liquid treatment system 

20 relates to water or wastewater treatment applications. Basic controllers, which are 
preferably feedforward and/or feedback controllers, can be combined into different 
control strategies. For example, the feedforward controller predicts the process 
changes and can optimize chemical consumptions. The fine tuning of chemical 
dosage may be done by the feedback controller. Controllers operate adaptively at 

25 large range of the process conditions. The change from one operating area to 
another is done by adaption model. The selection of proper control system depends 
on process requirements: a single feedback or feedforward controller can be 
sufficient; on the other end a cascade controller can be included for adapting the set 
point to special operating conditions. Preferably each control task can be realized as 

30 a single controller comprising a LE-controller and an adaptation model. 

Measurements can be direct or indirect real time or off-line measurements. Direct 
measurements are obtained by on-line instruments, which measure properties such 
as turbidity, suspended solids, conductivity, cationic demand etc. Indirect 
35 measurements (also called soft sensors) are calculated from on-line and/or off-line 
measurements e.g. by LE-model or regression model. Indirect measurements can 
replace expensive or missing on-line instrument measuring the properties of liquid. 
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5 

These measurements can be used in predictive mode indicating the change of 
operation area. 

In one embodiment the device arrangement for automatic dose control of chemicals 
5 in liquid treatment comprises one or more adaptation models and controllers, and an 
index describing the properties of liquid is used to modify the control surface of a 
linguistic equation (LE) controller, and the adaptation is defined with a linguistic 
equation model, to control the dosing of chemicals to the liquid by one or more 
controllers. 

10 

In another embodiment the device arrangement for automatic dose control of 
chemicals in liquid treatment comprises also an intelligent analyzer for monitoring 
the process. Intelligent analyzer is an implemented software module or device 
representing measurement handling routines. 

15 

In another embodiment the system comprises also means for data logging and 
means for connecting to data network (for example internet or LAN) for follow up 
of the controllers, monitoring, controlling and optimizing the process and updating 
the programs or means for connecting to telephone network. In still another 

20 embodiment the system is arranged to send an alarm message to a remote system, 
e.g. to a mobile phone, email, remote terminal or any other suitable system, when 
needed, for example in case of problem or for trouble shooting. In still another 
embodiment the adaptation system is based on remote operation, performed for 
example by using a data network. The parameters of adaptation models and 

25 controllers can be evaluated and updated when needed in a remote computer. New 
parameters are sent by data network to controller system which comprises for 
example a check routine for evaluating correctness of the parameters. 

The method according to the invention can also be used in dosing control of 
30 chemicals in retention, deposit control and foam prevention applications. Retention 
refers to the efficiency with which small particles (or additives) remain in the paper 
during its formation rather than staying with the white water. Deposits refer to the 
accumulations of material, coming from the water or suspended particles, onto 
wetted surfaces within a paper machine system or onto paper board. The method 
35 optimizes the efficiency of chemical treatment and the amount of chemicals. In 
addition, the wastes are minimized in the applications. As a result, the costs in the 
process are decreased. 
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The target for the chemical treatment efficiency itself depends on process 
requirements. Treatment can influence further on the productivity of the processes 
and quality of the product. Overdose of chemicals may produce an unacceptable 
product in the process. Target is to add the amount of chemical or chemicals that 
5 is/are needed to optimize the process performance and to avoid problems. 

The chemicals to be used in the method of the current invention can be any 
chemicals appropriate for the liquid treatment process used. Some examples of such 
chemicals are, but are not limited to, coagulants, flocculants, oxidants, reductants, 
10 adsorbents, dispersing agents, defoamers and biocides. 

Next the invention will be described in detail. The description will refer to the 
following drawings. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a schematic drawing of a water quality indicator for the case 
example. 

20 

Figure 2 is a graph depicting the purity index and scaled cationic demand. 

Figure 3 shows a schematic drawing of a feedback controller where coefficients Ci- 
C 3 are tuning parameters. RE mean real values and LE mean linguistic levels in the 
25 range of [-2,2]. 

Figure 4 shows a schematic drawing of an adaptive feedback controller based on 
quality of incoming water according to the soft sensor in water treatment 
application. 

30 

Figure 5 shows a schematic drawing of a generic adaptive model. 

Figure 6 shows a schematic drawing of a generic feedforward controller. 

35 Figure 7 shows a schematic drawing of a feedforward controller based on the 
quality of incoming water in water treatment application. 
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Figure 8 shows a schematic drawing of the dynamic simulator developed for the 
water purification unit. 

Figure 9a is a graph depicting the fluctuation of the water quality during the plant 
5 test on February 17. 2003 wherein the incoming water feed is constant. 

Figure 9b is a graph depicting the performance of the new LE feedback controller at 
the plant tested on February 17. 2003. 

Figure 9c is a graph depicting the simulated and the measured turbidities with test 
10 data of February 17. 2003. 

Figure 10a is a graph depicting the situation where the quality of the water is getting 
impure and the incoming water feed is constant. 

15 Figure 10b is a graph depicting the results from a test where the performance of 
current commercial feedback controller was tested on the plant. 

Figure 10c is a graph depicting the simulated and the measured turbidities with test 
data of March 10. 2003. 

20 

Figure lOd is a graph depicting the simulated performance of the new LE feedback 
controller on data of March 10. 2003. 

Figure 1 la is a graph depicting the situation where the quality of the water is getting 
25 more pure and the incoming water feed is constant during test on May 21-22. 2003. 

Figure 1 lb is a graph depicting the results from a test where the performance of 
commercial feedback and feedforward controllers were tested on May 21-22. 2003. 

30 Figure 11c is a graph depicting the measured and the simulated turbidities on the 
data of May 21-22. 2003. 

Figure lid is a graph depicting the simulated performance of the new LE feedback 
and feedforward controllers tested with the data of May 21-22. 2003. 

35 

Figure 12a is a graph depicting the situation where the amount of the treated water 
is halved during test on March 13. 2003, 
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Figure 12b is a graph depicting Hie results from a test where the performance of the 
new LE feedback and feedforward controllers was tested on March 13. 2003. 

Figure 12c is a graph depicting the measured and the simulated turbidities on the 
5 data of March 1 3 . 2003 . 

Figure 13a is a graph depicting the situation where the amount of the treated water 
is getting doubled and the quality of the water is rather constant. 

Figure 13b is a graph depicting the results from a test where the performance of 
1 0 current commercial feedback controller was tested on March 10-11. 2003 . 

Figure 13c is a graph depicting the measured and the simulated turbidities on the 
data of March 10-11. 2003. 

15 Figure 13d is a graph depicting the simulated performance of the new LE feedback 
controller tested with the data of March 10-11. 2003. 

Figure 14a is a graph depicting the situation where the amount of the treated water 
is getting halved. 

20 

Figure 14b is a graph depicting the results from a test where the performance of the 
current commercial feedforward and feedback controllers tested on June 4. 2003. 

Figure 14c is a graph depicting the measured and the simulated turbidities on the 
25 data of June 4. 2003 . 

Figure 14d is a graph depicting the simulated performance of the new LE feedback 
and feedforward controllers tested with the data of June 4. 2003. 

30 Figure 15a is a graph depicting the quality and the amount of incoming water on 
June 2-3. 2003. 

Figure 15b is a graph depicting the results from a test where the performance of the 
current commercial feedforward and feedback controllers were tested on June 2-3. 
35 2003. 



Figure 16a is a graph depicting the quality and the amount of incoming water on 
May 1-2. 2003. 
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Figure 16b is a graph depicting the results from a test where the performance of the 
current commercial feedback controller was tested on May 1-2. 2003. 

Figure 17a is a graph depicting the quality and the amount of mcorning water on 
5 February 5-6. 2003. 

Figure 17b is a graph depicting the results from a test where the performance of the 
current commercial feedback controller was tested on February 5-6. 2003. 

10 

DEFINITIONS 

As used herein the expression "the quality of liquid" describes the properties of the 
liquid which can be physical or chemical on-line measurements and/or on-line 

15 analysis or calculated values for example from one or more on-line/offline 
measurements and/or analysis of liquid, process conditions, consumption of 
chemicals, treatment result and target of the treatment efficiency. Examples of 
specific types of these kinds of calculated properties are quality index and purity 
index, which are explained below. The physical or chemical properties can be for 

20 example turbidity, conductivity, cationic demand, pH, amount of suspended solids, 
extractives or phosphorus. The quality of the liquid can be expressed as an index or 
factor, which can be in numeric or linguistic form. 

As used herein the term "quality index" (QI) describes the relative quality of the 
25 liquid to be treated by chemicals. Further the quality index describes the level of 
chemical(s) need to obtain the desired process performance. The quality index can 
be used in feedforward control. In addition, the quality index helps to detect the 
state of the process. 

30 As used herein the term "purity index" (PI) is a type of quality index, which 
describes the relative purity level of the incoming liquid to the liquid treatment unit. 
The purity index describes the need of chemicals to obtain the desired level of 
purity of the outgoing liquid. Current on-line measurements known in the art such 
as measurements of suspended solids do not well enough describe the level of 

35 chemical(s) need. The purity index correlates for example to the concentrations of 
soluble hemicelluloses and extractives substances. These concentrations can not be 
seen in on-line suspended solids measurement of the incoming liquid. 



WO 2005/022278 



PCT/FI2004/000501 



10 

As used herein the term "liquid" defines any substance containing sufficient amount 
of liquid phase to be used in a method of the current invention. The liquid can be for 
example substantially pure liquid or slurry which contains inorganic and/or organic 
particles or other substances, or something between these two options. Examples of 
5 liquids are, but are not limited to, water, waste water, industrial water, sludge or 
solids suspension, pulp suspension or any other liquid, such as a solvent or other 
chemical. 

As used herein the term "liquid treatment" means treatment of liquid by one or 
10 more chemicals, which are preferably dosed optimally to reach the efficiency 
required. The processes following the treatment process or the quality of the product 
or authorities give the target and limit for chemical treatment efficiency. Over dose 
or under dose of chemicals is avoided by proper control actions. Liquid treatment 
may contain several process stages for example coagulation, flocculation, oxidation, 
15 reduction, adsorption, absorption, foam prevention, deposit control and/or 
separations. 

Coagulation is the destabilization of colloidal particles brought about by the 
addition of a chemical reagent known as a coagulant. Fine particles in a suspension 
20 collide with each other and stick together. Usually the particles are brought near to 
each other by Brownian motion or by flow (Water Treatment Handbook, Vol 1 and 
Vol 2, 1991, Degremont). 

Coagulant is an inorganic (anion/cation) or organic (polyelectrolyte) chemical, 
25 which neutralizes the negative or positive surface charge (destabilization) of the 
impurities, such as colloidal particles. (Water Treatment Handbook, Vol 1 and Vol 
2, 1991, Degremont) 

Flocculation refers to the action of polymers in forming bridges between suspended 
30 particles or working by patch models, causing strong, relatively irreversible 
agglomeration. 

The reagent called a flocculant or a flocculant aid may promote the formation of the 
floe. Flocculant can be an inorganic polymer (such as activated silica), a natural 
35 polymer (starch, alginate) or synthetic polymers (Water Treatment Handbook, Vol 1 
and Vol 2, 1991, Degremont, Water Quality and Treatment, A Handbook of 
Community Water Supplies). 
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Oxidant can be any oxidising agent known in the art. Usually oxidant is a 
compound that spontaneously evolves oxygen, e.g. either at room temperature or 
under slight heating. Oxidants are generally used for bleaching and as part of 
biocides. Hydrogen peroxide and peracetic acid are examples. 
5 Adsorbent is a substance which has the ability to condense or hold molecules or 
ions coming out of aqueous solution and remaining on a surface. Bentonite is an 
example. 

Dispersing agent is a substance such as phosphate or acrylate which cause finely 
10 divided particles to come apart and remain separate from each other in suspension. 

Biocides are chemical additives designed to inhibit the growth of microorganism or 
to kill microorganisms, such as harmful slime-forming bacteria or fungi. Biocides 
are generally used in deposit control applications. 

15 

Defoamer (defoaming agent) is a substance used to reduce foaming due to gases, 
nitrogenous materials or proteins, which may interfere with processing. Examples 
are long chain fatty alcohols, organic phosphates, silicone fluid etc. 

20 As used herein "feedback controller" denotes any feedback controller known in the 
art. Generally the feedback control is an operation which, in the presence of 
disturbances, tends to reduce the difference between the output of a system and the 
set point (or an arbitrarily varied, desired state) by producing the actuating signal in 
such a way as to return the output to the desired value. A feedback control system 

25 is a system which tends to maintain a prescribed relationship between the output 
and the reference value by comparing these and using the difference as a means of 
control. (Stephanopoulos: Chemical Process Control. An Introduction to Theory and 
Practice, 1984, Prentice-Hall, p. 241-248) 

30 As used herein "feedforward" controller denotes any feedforward controller known 
in the art. Generally the feedforward control means control of effects of directly 
measurable disturbances or the process changes by approximately compensating for 
their impacts on the process output before they materialize. This is advantageous, 
especially in slow systems, because, in a usual feedback control system, the 

35 corrective action starts only after the output has been affected. (Stephanopoulos: 
Chemical Process Control. An Introduction to Theory and Practice, 1984, Prentice- 
Hall, p. 411-427). The feedforward controller can optimize also the ratio of 
chemical consumption when two or more chemicals are used. 
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A control system, which can adjust its parameters automatically in such a way as to 
compensate for variations in the characteristics of the process it controls, e.g. 
changes of desired steady-state operation, is called adaptive. In non-stationary 
processes the characteristics change with time. The optimal performance is 
5 maintained by adjusting parameters on the basis of an additional criterion. 
(Stephanopoulos: Chemical Process Control. An Introduction to Theory and 
Practice, 1984, Prentice-Hall, p. 431-438). In the systems of the current invention 
the adaptation is performed by the adaptation model. 

10 As used herein the term "adaptation model" defines the subsystem for modifying 
the control surface in changing operating conditions, which are defined for example 
with the quality index, e.g. purity index, and the process requirement, e.g. set point. 
As used herein the term "adaptation" includes also predefined modification 
procedures of the controller parameters, i.e. on-line self-tuning is not needed. 

15 Adaptation can be performed for example with linguistic equation (LE) models, 
which are especially beneficial in complex systems. The LE technique can combine 
data and expertise to represent non-linear interactions between variables as a model 
which is phenomenologically understandable. The adaptive LE controllers operate 
at large range of process conditions. The adaptation is based on detection of the 

20 process state. For simple systems the linguistic equation model can be replaced for 
example by fuzzy systems or conventional modeling techniques. 

As used herein "cascade" control system denotes any cascade control system known 
in the art. A cascade control system is a multiple-loop control system where a 

25 master controller controls the primary variable by adjusting the set point of a related 
secondary variable controller (slave controller). The secondary variable then affects 
the primary variable through the process. Disturbances arising with the secondary 
loop are corrected by the secondary controller before they can affect the value of the 
primary controller output. Usually there is only one master controller and only one 

30 slave controller, but some applications benefit from the use of more than one slave 
controller. (Stephanopoulos: Chemical Process Control. An Introduction to Theory 
and Practice, 1984, Prentice-Hall, p. 395-402). 

As used herein "intelligent analyzer" denotes an implemented software module or 
35 device representing measurement handling routines, e.g. indirect measurements, 
trend analyzers and detection of operating conditions. Intelligent analyzers can be 
used together with controllers. 
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DETAILED DESCRIPTION OF THE INVENTION 

The current invention is based on a non-linear linguistic equation (LE) controller 
whose control surface is modified adaptively by using the information obtained 

5 from the properties of liquid and the process requirements by means of a predefined 
adaptation model. The control system comprises one or more feedback and/or 
feedforward controllers for controlling the dosing of chemicals, adaptation models 
for performing the adaptation and optionally intelligent analyzers. A cascade 
controller may be used for special operating conditions where also process 

10 requirements need to be changed. 

The current invention can be applied in several applications of dosing of one or 
more chemicals in liquid systems such as, but not limited to, water, waste water and 
industrial water treatment, sludge treatment (dewatering, thickening), retention and 
15 deposit control. In these liquid systems chemicals are controlled by measurements 
that give response to chemical dosage. Liquid can be for example, but is not limited 
to, water, waste water, industrial water, sludge, suspension, liquid pulp suspension 
or any other liquid, such as any other solvent or chemical. 

20 

Linguistic equations 

The method used in the current invention is linguistic equation (LE) method which 
is a comprehensible tool for non-linear modeling and control applications. The 
25 linguistic equation (LE) contains of two parts: interactions are handled with linear 
equations, and non-linearity is taken into account by membership definitions. The 
general LE model can be presented as follows. 

AX + B = 0 , (1) 

30 

wherein the matrix X defines linguistic levels for variables and matrix A defines the 
direction and strength of the interactions between variables. Bias term B can be 
used to shift the model from the origin. Different operation points can be taken into 
account by adaptive features. The FuzzEqu Toolbox in the Matlab environment 
35 provides routines for developing the LE-models from the real data or based on the 
knowledge of the experts. This toolbox includes automatic generation of systems, 
model-based techniques and adaptation techniques (Juuso, 1999, Fuzzy control in 
Process Industry: The Linguistic equation approach. In: Verbruggen, H. B., 
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Zimmermann, H.-J., Babuska, R, editors, Fuzzy Algorithms for Control, 
International Series in intelligent Technologies, pp. 243-300. Kluwer, Boston, 
USA). 

Linguistic levels of the input variables are determined by means of membership 
5 definitions. These definitions scale the real values of variables to the linguistic 
levels with the range of [-2 - +2]. Membership definition comprises two 
monotonically increasing functions. The first function determines the linguistic 
values between -2 and 0, and the second function determines the linguistic values 
between 0 and 2. As an example these functions can consist of two second-order 

10 polynomials. Second order polynomials were used in the comparison examples. 
Membership definitions can be also defined by using the process expertise. 
(Jarvensivu etal, 2001 Intelligent control of a rotary kiln fired with producer gas 
generated from biomass. Engineering Applications of Artificial Intelligence, Vol. 
14, p. 629-653., Juuso, E.K. (2004). Integration of intelligent systems in 

15 development of smart adaptive systems. International Journal of Approximate 
Reasoning, 35, 307-337) 



20 



In the case of polynomial membership definitions, the linguistic level of the input 
variable j is calculated according to equation 2. 



2 with Xj > max(jc y ) 



2a) 



2a' t 



with Cj^Xj^ max(jc y ) 
with min(jc y ) < x y < c y 



(2) 



- 2 with Xj < imn(Xj) 



where aj,fcj,a) and bj are coefficients of the polynomials, Cj is real value 
corresponding to the linguistic value 0 and xj is the real value. Minimum and 
25 maximum values of the real data are min(x^ and max(^) corresponding to the 
linguistic values -2 and 2. After the linguistic level of the model output (X output ) is 
calculated according to equation 1, it is converted to real value of the output {x outpu ^ 
using the following equation: 



30 



a cvn«il X mumt + b ouaml^aiaui + C maui With A'„„„„, SO 



WO 2005/022278 



PCTYFI2004/000501 



15 

where a^ ut ,b; utpun a^ utput and b^ tput are coefficients of the polynomials and c output is 
the real value corresponding to the linguistic value 0. 



5 The example of water purification 

Next the current invention is described with examples where the method according 
to the invention is used in the purification of water in an industrial process. The 
following examples are intended to be exemplary and are provided to enlighten the 
10 invention and should not be considered as limiting. 

In one embodiment of the current invention a purification unit is connected in an 
internal water circulation in forest industry. The unit is used for removing of 
suspended solids and extractives from the process water in order to reuse the water 

15 in the process. The incoming process water contains generally a lot of suspended 
solids, extractives and anionic substances. Anionic substances are mainly soluble in 
water so they can not be removed by precipitation. However, these substances 
consume chemicals for neutralizing the cationic charge. The quality of incoming 
water to the flotation may fluctuate in a large range depending on the process 

20 conditions. The internal water circulation has also an effect on the quality of the 
incoming water. 

Two chemicals are fed to the incoming water before the basin. Coagulant 
neutralizes the charge on the surface of particles. It enables particles and colloids to 

25 coalesce and form little floes. Flocculant forms bridges between floes and increases 
the floe size. Part of the clarified effluent is recycled, pressurized, and saturated 
with air. The recycled water is fed to the flotation basin through the pressure release 
device. In the device the air is released from the water in the atmospheric pressure 
and micro bubbles are formed. In the reaction zone of the basin the air bubbles 

30 attach themselves to the floes. In the separation zone the bubble-floc-agglomerates 
float to the surface, where they are removed by mechanical scraping. 

Several water properties are measured continuously from the process. Concentration 
of suspended solids, conductivity, temperature and flow rate of incoming water are 
35 measured by on-line sensors. In addition, an on-line soft sensor estimates the 
proportional quality of the incoming water. Turbidity and pH are measured in the 
outgoing water. Measurements of suspended solids and turbidity are based on the 
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intensity of backscattering light with different wavelengths from the particles and 
colloids. 

The dynamic simulator was developed for the purification unit. Main purpose of 
5 the simulator was to create an environment for developing, testing and tuning the 
controllers and to create a tool for analyzing the operation of the process. Simulator 
contains a dynamic linguistic equation (LE) model for the flotation basin, 
controllers for two chemicals and a water quality indicator (Ainali L, Piironen M., 
Juuso E. (2002) Intelligent Water Quality Indicator for Chemical Water Treatment 
10 Unit. In Proceedings of SIMS 2002 - 43 rd Conference on Simulation and Modeling, 
Oulu, Finland, September 26-27). The simulator was built in Matlab Simulink. 

The faster effecting chemical, flocculant, is controlled by Pi-type feedback LE 
controller with an adaptation model. The change of control is calculated on the basis 
15 of the error and the change of error between the set point of turbidity and measured 
turbidity. More slowly affecting coagulant is controlled feedforward by a 
knowledge-based controller. Other control alternatives for coagulant are a steady 
state LE model and a constant value. 

20 

LE-model for the water purification unit 

Generally a LE-model for a water purification unit is generated from on-line data. 
Experimental design techniques are effective way to collect reliable on-line data. 
25 The modeling of the process is not necessary for implementation of controllers but 
it can ease tuning and testing of the control system in advance. 

In a case example central composite design was utilized in the development of the 
models and in order to evaluate the effect of chemicals dosage on the treatment 
30 results. Chemicals were tested in a large range to find real effects of the variables. 
Variables were tested at five different levels in order to evaluate the non-linear 
effects. The amount of experiments (9) is moderate compared to complete five-level 
factor design (25 experiments). 

35 The input variables of the model are the concentration of suspended solids (ss) 9 the 
previous value of the turbidity (turbiri) and the amount of chemicals (coag, floe). 
The output variable is the turbidity (turbout). The interaction coefficient of the 



WO 2005/022278 



PCT/FI2004/000501 



17 

variable j is ay Input variables are scaled to the linguistic values by equation 2. The 
linguistic level of the outlet turbidity (Y) is given by equation 4. 

Y = -lla turbout (a coag -X coag +a M *X floc +a„ X a -X turbUl ) (4) 

5 This basic model predicts the outlet turbidity quite well, when the quality of 
incoming water is typical and stable. Inaccuracy of the model increases when the 
quality of incoming water changes greatly e.g. due to the concentration changes of 
anionic substances at incoming water. These concentration changes cannot always 
be seen with current on-line measurements. 

10 

A dynamic multi-model was developed to avoid these problems occurring in the 
accuracy of the basic model. Multi-model comprises sub-models for different 
operation conditions. The structures of sub-models are same as those of the basic 
model. Same interaction matrix is included in all the sub-models. This is quite 

15 reasonable since the directions of interactions do not vary considerably between 
different operation points. The differences between the models are treated with 
membership definitions. The water quality indicator selects a proper sub-model for 
each operation condition. The multi-model is used when the controllers are 
developed and tuned in the simulator. The multi-model is not used in the calculation 

20 of a quality index for control. 

The quality index 

25 The quality index is a calculated value, which describes the quality of incoming 
water. Generally the quality of incoming water changes greatly in the treatment unit 
due to the concentration changes of impurities. These concentration changes cannot 
always be seen with current on-line measurements. In these situations the quality 
index is used for the monitoring of the water quality. 

30 

Generally a model based soft sensor calculates the quality index on the basis of on- 
line and/or off-line measurements and/or analysis, process conditions, the 
consumption of chemicals and/or purification result. The model can be any known 
model such as mechanistic model, regression model or LE model. Measurements 
35 and analysis such as physical or chemical properties of liquid, liquid and chemical 
flow rates are inputs of the model. Model results a new property describing the level 
of chemical need. 
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In the case example a water quality indicator calculates the purity index of incoming 
water. A basic model for the purification unit is a core of the water quality indicator. 
The indicator contains three stages: the basic model, the calculation of error 
between measured (on-line) and predicted outlet turbidity and the scaling of the 
error between -2 and 2 (Figure 1). If the error is positive the water quality is more 
impure than the average value. If the error is negative the water quality is more pure 
than the average. The average (normal) water quality was defined using on-line data 
for a long period (one month). 

The water quality indicator calculates the purity index of incoming water in the 
range of -2 to 2, which corresponds to the properties from extremely impure, 
impure, normal and pure to extremely pure. The cationic demand of the incoming 
water correlates strongly with the purity index of the incoming water. Figure 2 
shows that the cationic demand is decreasing and the incoming water is getting pure 
(the purity index increases from impure to extremely pure) at the same time. 
According to laboratory analysis the cationic demand characterizes principally the 
concentration of anionic hemicelluloses. The cationic demand was analyzed from 
samples by means of Mutec laboratory instrument. 

The feedforward controller 

The feedforward controller takes into account the process changes in advance. It can 
also optimize chemical consumptions. Input variables of the controller are for 
example direct and/or indirect calculated measurements, which describe the quality 
of incoming water and the need of the chemical dosage. 

The controller can be a static LE-model or a knowledge-based controller which is 
simplified form of the static LE- model. The static LE-model based controller can 
be expressed in the following equation, where Y is the linguistic value of the 
chemical dosage, a\ is an interaction coefficient and Xi is the linguistic value of 
input variables i, i=l . . .n. 

Y = (a x X, + a 2 X 2 + ... + a n X n ) (5) 

Scaling between real and linguistic values is carried out by membership definitions. 
The proper chemical dosages can be defined in the membership definitions. 
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If number of input variables is small and the effect of variables on output variables 
is linear the static LE-model can be simplified so that the scaling between real and 
linguistic values is not necessary. This kind of LE-controller is so-called knowledge 
based controller (Figure 6). It can be expressed in the following equation, where y is 
5 the value of a chemical dosage, y bas ic is a typical basic dosage, Xi is the value of an 
input variable and Wi is the weighting coefficient of an input variable i, i=l . . .n. 

ryw+^+^A (6) 

10 In a case example the slowly effecting chemical is dosed feedforward. The typical 
set point, x TSP , is achieved by the basic dosage when the incoming water is normal 
(purity index is 0). In the knowledge-based controller the differences in operating 
conditions are handled with the purity index of the incoming water, x Ph and the 
difference between the set point, x SPi and the typical set point Xjsp, (Figure 7). 



15 



y = y basic + Wpi x Pi + w S p(*sp - X KP) (7) 



The feedback LE-controller 

20 

The fine-tuning of chemical dosage is done by the feedback Pi-type LE-controller. 
The change of control is calculated on the basis of the error and the change of error 
between the set point of turbidity and measured turbidity (Figure 3). 

25 Input variables of the feedback LE-controller are a set point value and a 
measurement from outgoing water for example turbidity. This measurement should 
describe the treatment result of the purification unit. It should also give the response 
for the change of the chemical dosages. 

30 The LE-controller can be expressed in the following equations, where Y is a 
linguistic value, y is a real value and Cj are a tuning parameters i, i=l, 2, 3. Scaling 
between real and linguistic values is carried out by membership definitions. 

^change of control = C 3( C 1<% error C 2^ change of error) 

35 

y n&v control * previous control ^change of control 



(8) 
(9) 



WO 2005/022278 



PCT7FI2004/000501 



20 

The adaptive LE-controller 

An adaptive LE-controller contains the LE-controller and an adaptation model. The 
performance of the adaptation model is based on detection of the process state by an 
adaptation model. The adaptive model modifies the control surface according to this 
knowledge. The LE controller already handles non-linear operating conditions. 
Adaptation further extends the operating range of the controller without changing 
the controller. Therefore adaptive LE-controller can operate at large range of the 
process conditions. Each feedback and feedforward controller can have its own 
adaptation model. 

The adaptation model is based on the LE-model (Figure 5). The input variables are 
direct and/or indirect calculated measurements, which determine the process 
operation area. The equation calculates an adaptation coefficient, which is 
multiplied with the change of control. The adaption LE-model can be expressed in 
the following equation, where Y is the linguistic value of the adaptation coefficient, 
ai is an interaction coefficient of variables and Xi is the linguistic value of input 
variables i, i=l . . .n. 

Y = ~{a x X x + a 2 X 2 + ... + a n X n ) (10) 
Scaling between real and linguistic values is carried out by membership definitions. 

In a case example the adaptation coefficient (Y c ) is calculated by LE-model on the 
basis of the purity index (X Pind ) and the set point of turbidity (X^) (Figure 4). 

^cr 5 ~ ~ i a Pind X Pind + a sp X sp ) (11) 
a c 

The feedback controller with adaptation feature, which is used in the preferred 
embodiment of the invention, accommodates much better to the dynamic situations 
than a basic feedback LE controller. 

When the water quality is pure the adaptation coefficient decreases the change of 
control in order to prevent the oscillation of the turbidity. When the water quality is 
impure the adaptation coefficient increases the change of control in order to 
accelerate control actions. 
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Tuning methods 

Tuning of control parameters is generally done when the control system is 
implemented. It is not necessary to tune the controllers on-line because operational 
principle of LE-controllers is different than it is in self-tuning PDD controller, i.e. 
modification of the controller parameters is based on direct or indirect 
measurements (soft sensors, intelligent analyzers) in a predefined way, e.g. process 
models or heuristic knowledge. Tuning in advance facilitates fast operation in 
changing operating conditions. 

Tuning and testing of controllers is most easily and fast to do in the simulator in the 
case of complicated control systems. The simulator contains a dynamic process 
model and controllers. Generally experimental design is utilized in the development 
of the model, because models are often data based. 

Tuning and testing of controllers can also be carried out in the process if the control 
systems are simple for example only feedback LE controller without a predefined 
adaptive feature. 

In all these situations methods based on for example random search algorithms (for 
example chemotaxis algorithm, genetic algorithm) or grid search can be used 
effectively to tune the controllers in the process or in the simulator. Trial and error 
method can also be used to tune the controllers for example when the control 
system is simple. 

Cascade control 

Cascade controller systems may be needed to recover or prevent problems occurring 
in some special cases. Since the feedback LE controllers keep the controlled 
variables efficiently within acceptable limits, the set point should be fitted to the 
requirements of the other sub processes, e.g. properties of the treated liquid must be 
within specific limits. Cascade control can be also used for optimization of the 
treatment level in connection to the overall process. The cascade improves both 
feedback and feedforward control. 
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EXAMPLES 

Comparison test examples 

5 

Comparison tests between commercial self-tuning PID-controllers and the new 
adaptive LE-controllers of the current invention were carried out during spring 2003 
in a water purification unit The unit was connected to an internal water circulation 
in forest industry. The incoming process water contains suspended solids, 

10 hemicelluloses, extractives and inorganic substances. The quality of incoming water 
to the purification unit may fluctuate in the large range depending on the process 
conditions. The internal water circulation has also an effect on the quality of the 
incoming water. The unit is used for removing of suspended solids and colloidal 
substances from the process water. In the following tests the concentration of 

15 suspended solids in the incoming water, the flow rate of incoming water and water 
turbidity in the outgoing were measured by on-line sensors. Measurements of 
suspended solids and turbidity are based on the intensity of backscattering light with 
different wavelengths (600 and 900 nm) from the particles and colloids. 

20 Two polyelectrolytes (Cheml, Chem2) were fed to the incoming water before 
separation stage. Cheml, which is a coagulant, neutralizes the charge on the surface 
of particles. It enables particles and colloids to coalesce and form little floes. 
Chem2, which is a flocculant, forms bridges between floes and increases the floe 
size. Floes are removed from the water in separation stage. Purification result is 

25 followed by on-line sensor of outgoing turbidity. The chemicals have been dosed 
either by self-tuning PID controllers or the new adaptive LE-controllers according 
to the current invention on the plant at a time. 

The basis of the commercial controller used in a comparison is a self-tuning PID 
30 controller designed exclusively for polymer control. Primary control variable is 
usually measured with on-line turbidity or suspended solids sensors. The 
performance of the controller is occasionally limited in rapid variations of operation 
conditions (e.g. fast turbidity changes) as the adaptation mechanism cannot work 
properly. 



35 



The commercial controllers and the new controllers of the current invention were 
tested separately on the plant. Exactly similar operation conditions for controllers 
cannot be repeated because of complex internal water circulation. Several sub- 
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processes have also a great effect on the amount and quality of water. Therefore 
process simulations with new controllers have also been used for the performance 
comparison. The simulator developed for the purification unit contains a dynamic 
linguistic equation (LE) model for the purification unit, two controllers for dosing 

5 of chemicals and a soft sensor describing the relative quality of incoming water 
(Figure 8). The model calculates outlet turbidity on the basis of chemical dosages, 
suspended solids at incoming water and flow rate of incoming water. The dynamic 
model has been tested with on-line data collected from the plant. According to these 
tests, the generated dynamic models provided a good prediction for the outlet 

10 turbidity (Figures 9c, 10c, 11c, 12c, 13c and 14c). 

The water quality indicator is a model based soft sensor. It calculates the purity 
index of incoming water in the range of -2 to 2, which corresponds to the properties 
of water from extremely impure, impure, normal, and pure to extremely pure. 

15 Typically the quality and the amount of incoming water can change a lot in a broad 
range. Quite often the changes are very fast. Keeping the purification result stable in 
the dynamic situations is very challenging for controllers. In the following examples 
the performance of the commercial and the new adaptive LE controllers of the 
present invention are compared to each other in changing operating conditions. 

20 

Definitions for parameters in Figures of comparison tests 

"Incoming Flow" (%) describes the relative amount of incoming water. The value 
25 of 100 % corresponds to the maximum capacity of the purification unit. 

"Purity Index" describes the quality of incoming water at the range -2 and 2 (from 
extremely impure to extremely pure water). 

30 'Turbidity" (NTU) describes the amount of colloidal substances in outgoing water. 

"Set Point" (NTU) describes the desired outgoing turbidity. 

"Scaled Dosages" describe relative amount of chemicals. The value 1 corresponds 
35 to the maximum amount of the chemicals. 
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Example 1. Change of incoming water quality 

Example la. Performance of the new adaptive feedback controller according to 
5 the invention when the water quality is changing. 

New adaptive feedback LE-controller can react to the change of water quality fast 
and can keep the set point value. There is not offset between the measured and set 
point turbidity. The change of the incoming water quality has an immediate effect 
10 on the chemical dosage because the feedback controller reacts fast to the error 
between turbidity and set point. 

In the example the purity index varies between 1 and 1.6 (Figure 9a). When the 
water is getting first slightly more pure the chemical dosage is decreasing 
15 immediately (purity index from 1.2 to 1.6). After that the impurity of water starts to 
increase slightly (purity index from 1.6 to 1.0). The feedback controller increases 
chemical and can keep the set point value (Figure 9b). 

20 Example lb. Performance of the commercial controller when the water quality 
is changing. 

The commercial feedback controller cannot always keep the set point value when 
the quality of incoming water starts to change, especially in fast changes. Controller 
25 starts to increase or decrease the chemical feed too slowly. As a result the error 
between measured and set point turbidity can be significant. Often a constant 
deviation from set point value can also remain. 

In the example the incoming water is very pure for 5 hours (purity index is 2, very 
30 pure) (Figure 10a). Then the quality of water is getting impure and approaches the 
normal water quality level (the purity index decreases from 2 to 0). Figure 10b 
shows that the commercial feedback controller keeps the set point when the water 
quality is stable. During the change the controller increases chemical dosage too 
slowly. As result the chemical dosage is not enough to keep the outlet turbidity in 
35 the set point value. 

According to the simulation results the new feedback controller can react to the 
change of water quality much faster and can keep the set point value (no offset 
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between the measured and set point turbidity, Figure lOd). Figure 10c shows that 
the dynamic model can predict the outlet turbidity well during this trial, 

5 Example lc. Performance of commercial feedforward and feedback 
controllers when the water quality is changing. 

The commercial feedforward and feedback controllers can keep the set point well 
when the quality of the incoming water is stable. Turbidity can deviate from the set 
10 point quite a lot when the water quality changes as it can be seen in Figures 11a and 
1 lb. According to the simulation results the new controllers can approach to the set 
point smoothly in the same situations (Figure lid). Figure 11c shows that the 
simulation model is reliable and predicts the measured turbidity well. 

15 In the beginning the quality of the incoming water is normal (purity index is 0). 
Small peaks in purity index (e.g. from 0.5 to 0, from 0.5 to -0.25, from 0.8 to 0.1) 
lead to peaks in turbidity with commercial controllers, but not with the new 
controllers. When incoming water is getting pure (purity index from -0.25 to 1) the 
commercial controllers can not react to changes in turbidity fast enough. The new 

20 controllers work better in this situation according to the simulations. The new 
controllers approach the set point faster and the deviation from set point is smaller. 
In addition, the commercial controllers lead to offset from set point during small 
long-run decrease of impurity level at the end of test (purity index decreases from 1 
at 8 a.m. to 0.7 at 6 p.m.). This offset is not observed with new controllers. 

25 

Example 2. The amount of water is changing 

30 Example 2a. Performance of both new controllers when the amount of water 

is changing (half of water feed) 

The new controllers can react to the halving of incoming flow quite well (Figure 
12a). Controllers can keep the set point value and there is no offset between the 
35 measured and set point turbidity as it can be seen in Figure 12b. 



In the example the set point of the turbidity is decreasing at the same time as the 
incoming flow is halving because of the process demands (optimization of sludge 
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flow). Before the change the controllers keep the set point very well. During the 
change of the flow rate the dosage of Cheml increases temporarily quite high 
because the chemical pump reacts slowly to the change of water flow rate. As result 
the outlet turbidity decreases also temporarily. Otherwise the outlet turbidity 
5 follows well the set point. 

Example 2b. Performance of commercial feedback controller when the 
amount of water is changing (double water feed) 

10 

The change of the water feed disturbs the operation of the purification unit. 
Doubling of the water feed is more demanding for the controllers than halving of 
the water feed as the amount of water increases in few minutes. This leads to a rapid 
peak in outlet turbidity due to large amount of water to be purified. Neither 
1 5 commercial nor new controllers can fix this peak completely. 

In the example the commercial feedback controller doses the chemical (Chem2) 
during the doubling of the water feed (Figure 13a). During the change Ihe feedback 
controller reacts too weakly to the peak of turbidity. As result the turbidity increases 
20 quite high (475 NTU). Figure 13b shows that the commercial feedback controller 
adds the chemical dosage too slowly after the change so the offset between set point 
and measured turbidity remains. 

According to the simulations the new feedback controller is working better in this 
25 situation (Figure 13d). The controller reacts more strongly to the peak of turbidity. 
As a result the turbidity increases only to the value of 350 NTU during the change. 
After the change the feedback controller keeps the set point value very well. 

30 Example 2c* Performance of both commercial controllers when the amount 

of water is changing (half of water feed) 

Halving of the water feed does not disturb the operation of the purification unit as 
much as doubling of the water feed. So this situation should be easier for the 
35 controllers. When both commercial feedback and feedforward controllers dose the 
chemicals during the change of water feed the changes in turbidity are fixed by 
feedforward controller. As result the turbidity can deviate from the set point values 
quite a lot. 
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In the example both commercial controllers dose the chemicals. During the halving 
of the water feed the feedforward controller reacts to the decrease in turbidity 
(Figure 14a and 14b). Cheml effects on turbidity more slowly than Chem2. That is 
why the turbidity decreases below the set point. After the change the controllers 
5 keep the set point momentarily. When the water is getting impure the controllers do 
not react to the increase in turbidity. As a result the offset between turbidity and set 
point value remains. 

The new feedback controller reacts to the change in turbidity (Figure 14d). As result 
turbidity does not decrease as low as with commercial controllers. The feedforward 
10 controller reacts to the change of water quality. 



Example 3. Disturbances in the use of the operation of commercial controllers 

15 Disturbances can sometimes appear in performance of commercial feedforward and 
feedback controllers, especially when the purification unit is starting up or the 
incoming flow rate is changing. In these examples the dosages of the chemicals can 
remain at the low or high limit. 



20 

Example 3a, The flow rate of Chemical 1 remains at the low limit 



The dosage of Cheml remains at the low limit during the start-up of the purification 
unit (Figure 15b). The feedback controller can keep the set point when the incoming 
25 flow rate of water was low (approximately 40%). When the incoming flow rate was 
doubling, the water quality was getting more impure (Figure 15a). The feedback 
controller can not keep alone the turbidity at the set point. Turbidity was increasing 
above the set point because of the low dosage of Cheml. 



30 

Example 3b. The flow rate of Chemical 2 remains at the high limit. 

The dosage of Chem2 remains at the high limit in the start-up of purification unit 
(Figure 16a and 16b). As a result the feedback controller can not decrease the 
35 dosage. Therefore outlet turbidity was below the set point. The disturbance was 
fixed manually by decreasing the dosage of Chem2. 
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Example 3c. The flow of Chemical 2 remains at the low limit. 

Figure 17b shows that the dosage of Chem2 remains at the low limit when the 
amount of incoming water is halving in few minutes (Figure 17a). Temporarily the 
5 turbidity of water is decreasing below the set point. The controller decreases the 
amount of chemical, but for some reason the dosage remains at low limit for several 
hours although the turbidity increases even over the set point. The amount of 
Chem2 is so small that the purification efficiency is declining and the turbidity is 
increasing above the set point. However, the controller does not increase the amount 
10 of Chem2. This disturbance of the feedback controller was fixed manually by 
increasing the dosage of Chem2 temporarily. 

This invention has been described with an emphasis upon some of the preferred 
embodiments and applications. However, it will be apparent for those skilled in the 
art that variations in the preferred embodiments can be prepared and used and that 
15 the invention can be practiced otherwise than as specifically described herein within 
the scope of the following claims. 



